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^ ' Abstract 

Q_i, We present an informal review of results on asymptotics of orthogonal polynomials, 

stressing their spectral aspects and similarity in two cases considered. They are 
polynomials orthonormal on a finite union of disjoint intervals with respect to the 
Szego weight and polynomials orthonormal on R with respect to varying weights and 
having the same union of intervals as the set of oscillations of asymptotics. In both 
cases we construct double infinite Jacobi matrices with generically quasiperiodic 
coefficients and show that each of them is an isospectral deformation of another. 
! Related results on asymptotic eigenvalue distribution of a class of random matrices 

of large size are also shortly discussed. 

o 
m . 

1 Introduction 

The goal of the paper is to discuss a link between asymptotics of a class of orthogonal 
J ■ polynomials, in particular polynomials with respect to varying weights (see e.g. |39j). 

and the Jacobi matrices with quasi-periodic coefficients, seen mostly as a particular case 
of ergodic finite-difference operators. The theory of this class of operators owes a lot to 
B. Simon, starting from an early book ^Hj till just appeared impressive |3b| . The link 
became clear while the author was reflecting on applications of the asymptotic formulas, 
found in the remarkable paper by P. Deift et al |15j . to certain problems on the eigenvalue 
distribution of a class of random matrices, known as unitary invariant matrix models. 
This is why we would like to begin from a discussion of random matrices, despite that 
the link can be described without a recourse to random matrices. 
Consider n x n Hermitian random matrices 



M n = {M jk e C, M kj = M jk }l k=1 , (1.1) 

whose probability law is 

P(dM n ) = Z- 1 exp{-nTrV(M n )}dM n . (1.2) 

Here Z n is the normalization constant, V : M. — > M. + is a continuous, bounded below 
and growing at infinity function (think about a polynomial of an even degree, positive at 
infinity, see also (|2.11jl ). and 

n n 

dM n = Y[dM jj [ [ d^tM jk d^M jk . (1.3) 

j=l ^<3<k<n 
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This class of random matrices arises in a number of fields of mathematics and physics 
(see e.g. reviews [13 EH [23 123 123 1211 and references therein). A considerable amount of 
corresponding problems can be described in terms of the Normalized Counting Measure 
of eigenvalues (NCM), defined as the relative to n number of eigenvalues of M n , falling 
into a given set AcR: 

N n (A) = #{A[ n) g A, Z = l,...,n}/n, (1.4) 

where 

{A[ n) }r=i (i.5) 

are eigenvalues of M n . 

We note that similar measures arise in spectral theory of ergodic operators, and B. 
Simon did a lot of excellent work on the measures, whose limit as n — > oo is known there 
as the Integrated Density of States. 

It will be convenient to consider a bit more general object, known as a linear eigenvalue 
statistics and defined via a test function ip : R — > R and eigenvalues (jl.5|) : 

n 

N n [^]=n- l Y,^i l) )- (1.6) 
i=i 

We obtain (jl.4|) by setting ip = xa, where xa is the indicator of A. 

Here are three basic quantities of the Random Matrix Theory, related to eigenvalue 
statistics and widely studied, especially as n — > oo. 

(i) Expectation of N n [ip] with respect to (jl.2|) - (|1.3j) : 

N n [(f] = E{N n [ip}}. (1.7) 



(ii) Covariance of N n for two test functions (pi.2- 

Cov^^i],^^]} (1.8) 
= EiNMNM} - E{JV n [^i]}E{JV n [v?2]}. 

(iii) Gap probability 

E n (A)=P{N n (A) = 0}. (1.9) 



It is a remarkable observation by Gaudin, Mehta, and Dyson of the early 60th (see e.g. 
[23] ) that the above quantities can be expressed via the orthonormal polynomials {p\ }i>o 
with respect to the weight 

w n (X) = e- nV(x \ (1.10) 



w n (\)p\ n) p£>d\ = 6 lmj l,m = 0, 1, ... (1.11) 

Here and below integrals without limits denote integrals over R. We will call V potential. 
Denote 

ti n \\)=wW{\)pi n \\), (1.12) 
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and 

n-1 

K n {X^) = Y,^ n) W4 n) ^)- (1-13) 

is called the reproducing kernel of the orthonormal system {^ n ^(A)}j>o, and (K n (X, A)) -1 
is known in the approximation theory as the Christoffel function |33J EH] • 

To distinguish these polynomials from the traditional ones for which the weight does 
not contain the large parameter n, the polynomials (|1.1U|) - (jl.l2|) are called the orthogonal 
polynomials with respect to varying weights (see e.g. [3T?]). 

We have for JTZJ) - lO|) : 

N n [<p] = j <p(\)p n (\)d\, Pn (X) = n- l K n {\X), (1.14) 

see p5] . 

Cov{NM, N n {<f 2 ]} = ^ j j A ^ A ^ K 2 n (X 1 ,X 2 )dX 1 dX 2 , (1.15) 



where 

Aip _ ip(Xi) - ip{X 2 ) 
AX~ Ai - A 2 

see (BO], an d 



;i.i6) 



E n (A) = det(l - K n (A)), (1.17) 
where K n (A) is the integral operator 



(K n (A)f)(X)= [ K n (X,p)f(p)dp, Xe A, 



(1.18) 



see |25! • 

We will need two more basic facts on orthonormal polynomials. The first is the Jacobi 
matrix J^ n \ associated with the polynomials via the r.h.s. of the three term recurrence 
relation: 

ApW(A) = r^pffiCA) + st ] p\ n \X) + r&p&(A), r_ t = 0, (1.19) 

i.e., 

7(n) _ r tWioo j(n) _ (n)p , (nU , (n) p /-, 9f| N 

The second is the Christoffel-Darboux formula 

A — /i 

Formulas (fTTT^I) - ffTTTH|) and (ITT^Tjt show that the asymptotic form of (JU7J) - (JH§ as 

n — > oo is determined by that of V'n-i an d These were found by Deift et al ^3] i n 
the case, where V is real analytic (see also [SI IS] ) - In the same paper the limit of (|1.17j) 
was found in the so-called local asymptotic regime (see e.g. [27| for its definition, and [3U] 
for another derivation of this result). 

The author of this paper has applied the results of [H^ to find an asymptotic form 
of the covariance (|1.8|) . and to prove an analog of the central limit theorem for linear 
eigenvalue statistics (jl.6|) . It turned out that this requires a bit more information on the 
asymptotic formulas of [15J, and leads to certain objects, related to spectral theory of 



3 



quasiper iodic Jacobi matrices. This is discussed in the paper. Asymptotics of covariance 
and an analog of the central limit theorem for linear statistics of eigenvalues of random 
matrices (jl.2|) - (jl.3J) will be published elsewhere |2~Sj . 

The paper is organized as follows. In the next section asymptotic formulas for the 
"ordinary" polynomials orthogonal with respect to weights whose support is a union of 
q > 1 disjoint intervals are shortly discussed following papers (21 ESI HOj • We present 
then asymptotics found in ^S] for polynomials, orthogonal with respect to varying weights 
for the case, where their oscillatory part is the same union of q intervals. In Section 3 
we introduce quasiperiodic Jacobi matrices, associated with the both asymptotics and 
discuss links between the matrices. We argue that they are so-called "finite-band" Jacobi 
matrices, widely known in the theory of integrable systems, that they are related by an 
isospectral deformation, and consider a particular case of polynomial potentials in ()1.10p . 
where corresponding Jacobi matrices are periodic. In Section 4 we present a collection 
of facts on asymptotic eigenvalue distributions of random matrices, that can be written 
in the terms of the above Jacobi matrices. In Appendix we give a direct proof of the 
isospectrality of the Jacobi matrices related to asymptotics of both classes of orthogonal 
polynomials. 

2 Asymptotics of orthogonal polynomials 

2.1 Ordinary orthogonal polynomials 

Consider first the case, where the weight w does not depend on n. In this case all the 
quantities, related to orthonormal polynomials, do not depend on the super-index (n). 
We will denote the polynomials and related quantities by the same symbols as in (jl.llj) 
- (|1.13|) and (|1.19|) - ()1.21|) but without the super-index (n). Assume that the support a 
of the weight is a finite union of disjoint finite intervals: 
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a = 

i=i 



\^J[ai,bi], — oo < a\ < b\ < ... < a q < b q < oo. (2.1) 



Denote .Mi(cr) the set of non-negative unit measures on a and consider the quadratic 
functional 

£er[m] = — / log |A — fi\m(d\)m(djj), m G M\{a). (2.2) 

JaXcr 

The functional possesses a unique minimizer v (the equilibrium measure for a): 

min £ a [m) = £ a [u\. 

This is a standard variational problem of potential theory, that admits a simple electro- 
static interpretation in which m is a distribution of positive charges on a conductor a and 
v is the equilibrium distribution of charges. 

It is known (see e.g. jSl]) that the problem is equivalent to the relations 

-2 / log | A - n\u(d/i) = -l a , Ago-, (2.3) 

J a 

-2 / log |A - fi\v(dn) > -l a , A G R\ a, (2.4) 

J a 
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that are the Euler-Lagrange equations for (|2.2|) . The quantity —l a /2 is known as the 
Robin constant, and e 1 "^ 2 is the logarithmic capacity of a. 
Set 

z/(A) = K(A,oo)), (2.5) 

and 

a = a i = v{ai+i). (2.6) 

With this notation the asymptotic formulas of [Zj, Theorem 1, (see also [210123110]) for 
analogs of orthonormalized functions (|1.12j) . corresponding to a n-independent weights w, 
can be written as follows. 

Assume that the weight satisfies the Szego condition 

logw £ L x {o, u). (2.7) 

Then there exist the functions V a : a x T^ 1 -> R+, and Q a : a x T^ 1 -> R such that if 
A belongs to the interior of a, we have 

^ n (A) = (2D CT (A,na)) 1 / 2 (2.8) 
x cos (imi>(\) + <? CT (A, na) j + o(l), n — > oo, 

where the remainder vanishes in the L 2 (cr)-norm, and 

na = (na%, ...na q -i). (2.9) 

Besides, there exist functions TZ a : T 9-1 — > E + , and <S CT : T 9_1 — > R, such that the 
coefficients {77, s;};> of the corresponding Jacobi matrix J a (|1.2(J|) . that does not depend 
on the super-index n in this case, have the following asymptotic form 

r n = 1Z a (na) + o(l), s n = S a (na) + o(l), n — > 00, (2.10) 

Note that to find (|2.10|) one needs weaker asymptotics of p n (A), those for A outside a. 

Functions V a , Q a , lZ a , and S a can be expressed via the (q — 1) - dimensional Riemann 
theta-function (see e.g. formula (jA.8|) below), associated with the two-sheeted Riemann 
surface. The surface is obtained by gluing together two copies of the complex plane slit 
along the gaps (61, a 2 ), (b q -i, a q ), (b q , a\) of the support of the measure u, the last gap 
goes through the infinity [Jj. For another form of V a , Q a) lZ a , and S a see jSSl ESI- 

The case q = 1 in ()2.1|) of polynomials orthogonal on a single interval dates back to 
S. Bernstein, Szego and Akhiezer [37]. 

The components of the vector a = {ai}1~l are rationally independent generically in a, 
thus the sequences {V a (X, n«)} ng z, and { Q a ( A , na)} n£ % for any fixed A and the sequences 
{7Za(na)} ne z, and {S a {na)} n( zi are quasiperiodic in n (see [JJ [TTJ [22 E3)- As an early 
precursor of this fact we mention a result by Akhiezer pQ, according to which if a consists 
of two intervals, then a certain characteristic of corresponding extremal polynomials of 
degree n can be expressed via the Jacobi elliptic functions as n —>■ 00. As a result the 
characteristic does not converge as n — > 00 but has a set of limit points that fill a specific 
interval generically in the intervals lengths in ()2.1|) . 
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2.2 Orthogonal polynomials with respect to varying weights 

Let V : R — > R + be real analytic and such that 

lim V(X)/\og(X 2 + 1) = oo. (2.11) 

| A|— >oo 

Consider orthonormal polynomials (jl.lOj) - ()1.12j) . To describe their asymptotics we 
introduce the functional (cf ()2.2|) ): 

£ v [ m ] = -J J log |A - fi\m(d\)m(dfj,) + J V(X)m(dX), (2.12) 

defined on the set A4i(R) of non-negative unit measures on R. 
The functional 1)2.12)1 possesses a unique minimizer N 

mjn £ v [m]=£ v [N\- ( 2 -13) 

The variational problem, defined by ()2.12jl . goes back to Gauss and is called the minimum 
energy problem in the external field V (see recent book [31] for a rather complete account 
of results and references concerning the problem). The unit measure N minimizing ()2.12|) 
is called the equilibrium measure in the external field V because of its evident electrostatic 
interpretation as the equilibrium distribution of linear charges on the ideal conductor 
occupying the axis R and confined by the external electric field of potential V. We stress 
that the corresponding variational problem determines both the (compact) support a of 
the measure and the form of the measure. 
The problem is equivalent to the relations 



$(A) = -l v , X E a, (2.14) 



$(A) > -l v , X e R\a, (2.15) 

where 

$(A) = V(A) - 2 ! log |A - n\N{dn). (2.16) 

J a 

This should be compared with the variational problem ()2.2|) of potential theory, where 
the external field is absent but the support o is given. This problem can be viewed 
as a particular case of ()2.12|) . corresponding to a sequence of potentials approaching 
V = Xo- 1 ~ 1) where Xa is the indicator of a. 

The minimum energy problem in the external field arises in various domains of analysis 
and its applications ^JOUEni- We will use here a link with Random Matrix Theory. It 
was argued by Wigner in the 50th (see j2H] for references and discussions) , and is shown in 
[TIH EI] that the measure N n of ()1.7)) converges weakly as n — > oo to the unique minimizer 
N in (|2.13|) . Moreover, the random measure ()1.4j) converges weakly to iV with probability 
1 as n — > oo. 

Assume that V is such that the support of is a union of q disjoint intervals as in 
f!2.1|) . Introduce the non-increasing function (cf ()2.5|) ) 

JV(A) = JV((A,oo)), (2.17) 
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and the (q — 1) - dimensional vector (cf ([2.6(0 

/3 = {A}?=i\ fl=JV(aj+i). 



(2.18) 



With this notation the asymptotics for orthogonal polynomials with varying weight found 
in [15] . Theorem 1.1, can be written as follows. There exist continuous functions IV : 
a x T 9_1 -»• R+, and Q v : a x T 9 " 1 -> K, and < r < 1 such that if A belongs to the 
interior of the support a ()2.1|) of N, we have (cf ([2.8(0 : 

#)(A) = (2P y (A,n/3)) 1 / 2 (2.19) 
x cos h7miV(A) + Qy{\ n0) \ + 0(n _T ), n -> oo, 

where ^ (A) is defined in f|1.12|) . and n(3 = (nfii, ...,n/3 g _i). If A belongs to the exterior 
of a, then ?/4 n ^ decays exponentially in n as n — > oo. 

Similar asymptotic formulas are valid for coefficients of the Jacobi matrix of (|1.2(J|) . 
Namely, according to ^5J, formula (1.64), there exist continuous functions IZy '■ T 3-1 — > 
R+ and S v : T^ 1 -»• R such that we have (cf flTTDD ) 

r W = fty(n/3) + 0(n- T ), sjj»> = <Sy(n/3) + 0(n- T ), n -> oo. (2.20) 

It will be argued below that the functions Vv,Qv,T^v, and Sv differ from the functions 
T^a-i Ga, 7£<t and S a of formulas ()2.9|) and ()2.10j) of the previous subsection only by a shift in 
the argument. Hence the main difference in asymptotic formulas of the previous and this 
subsection is that in the former the "rotation number" and the frequencies are determined 
by the measure v (see (12 .5j) - ([2. 6(0 . minimizing the functional (J2.2)) . while in the latter 
these quantities (see ()2.17|) - ()2.18[) ) are determined by the measure N, minimizing the 
functional (l2~T2l . 



3 Quasiperiodic Jacobi matrices 
3.1 Ordinary orthogonal polynomials 

Consider orthonormal polynomials with respect to the weight, whose support is a union 
of q disjoint intervals (J2.1j) . Denote J a the semi- infinite Jacobi matrix, associated with 
the polynomials and let {r h si}i> be the non-zero coefficients of J a (see ([1.19(1 - ([1.20(1 . in 
which the super-index (n) is omitted). Introduce the double-infinite Jacobi matrix J CTi7l , 
setting 

_ ( r k+n , k>-n, _ J s k+n , k>-n, , . 

r ^-\ o, k<-n, s ^-\ o, k<-n. 

We will denote by the same symbol J a , n the selfadjoint operator in / 2 (Z), defined by the 
matrix. 

Now the asymptotics (|2.10() allow us to define a family of "limiting" double- infinite 
Jacobi matrices and corresponding selfadjoint operators in Z 2 (Z). Assume for the sake of 
definiteness that the components of the (q — l)-dimensional vector a of ([2.6(1 are rationally 
independent. Then for any x — (x%, x q -i) G T 9 ^ 1 there exists a subsequence {rii(x)}i>i, 
such that 

lxm{rii(x)ai} = x h I = 1, ...,q- 1, (3.2) 
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where {t} denotes the fractional part of t G R. This and (|2.1U|) imply that for any k G Z 
we have 

lim r n j x )+k = TZ a (ka + x), lim s n j x ) +k = S a (ka + x). (3.3) 

i— >oo i— too 

In other words the sequence {J a ,m(x)}i>i of selfadjoint operators, defined in / 2 (Z) by the 
double infinite Jacobi matrices with coefficients (jH.ljl . converges strongly to the operator 
in Z 2 (Z), defined by the double - infinite Jacobi matrix J a (x) with coefficients 

n a (ka + x), S a (ka + x), k <E Z. (3.4) 

The matrices J CT (x) arise in spectral theory and integrable systems [T71 1221 EH] and is 
known there as finite band Jacobi matrices. 

Write the three-term recursion relation for J CT n : 

r n+k ip n+k+1 + s n+k ip n+k + r n + k -Y*\)n+k-\ = Xip n+k , k > —n, A G o. 

Setting here n = rii(x), using asymptotics ()2.8jl . and taking into account that in the 

obtained asymptotic equality the coefficients in front of "fast oscillating" expressions 

cos(7rn_ ? -(x)z>'(A)) and sm(7rrij(x)u(X)) should be both zero at the limit i — > oo, we find 
that for any A, belonging to the interior of a, the sequences 

{{V a {\, ka + x)) 1/2 cos(vrz/(A)A; + G a (ka + x))} k& , (3.5) 

and 

{{V c {\ ka + x)) l/2 sin(7rz/(A)A; + Q a {ka + x))} k& (3.6) 

satisfy the limiting three term recurrence relations, defined by the coefficients (|3.4|) . In 
other words, sequences ()3.5j) and ()3.6|) are generalized eigenfunctions of J ff (x) for every 
A, belonging to the interior of a. 

Note now that by general principles [3] the resolution of identity £j CT of the initial 
Jacobi matrix J a is 

{£j„(d\)) ik = ip j (X)ip k (X}dX, j, k>0, (3.7) 

in particular 

f (£ J<7 (d\)) jk = 5 jk . (3.8) 

Hence, the resolution of identity £j of J^, defined by ()3.1|) . is 

(£ (d\)) = I ^n+j(X)tfjn+k(X)dX, j,k>-n 
1 Ja ' n{ >hk \ 0, otherwise. 

This and asymptotics f)2.8j) yield for the weak limit of the above projection-valued measure, 
the resolution of identity £j a ( x ) of J a {x): 

{£ Mx) {dX)) jk = Xa(X)(V a (X,ja + x)V a (X, ka + x)) 1 / 2 (3.9) 
x cos [irv(X)(j — k) + Ga-{joi + x) — Q a {ka + x)jdX, j, k G Z, 

where Xa is the indicator of a. Denoting ifj(x) = 7rz/(A)j + Q(ja + x), j G Z, we can write 
the cosine above as cos<pj(x) cosipk(x) + sin^^x) sm<fk(x)- This shows that the r.h.s. of 
()3.9|) is the linear combination of ()3.5|) - (J3.6|) . Besides, the equality 



/ {.£j a (x)(dX)) jk = 5 jk , j, k G Z 

J a 



that can also be proved by the limiting transition rii(x) — > oo in (|3.8|) . implies that the 
union of the sequences (|3.5|) - (|3.6|) for all A of the interior of a forms a complete system 
in 1 2 {Z). 

Introducing 

* j (\,x) = e im,( ' x)i u j (\,x), (3.10) 

where 

Uj (X,x) = U(X,ja + x), U{\x) =V 1 J 2 (\,x)e iG ° {x > x \ (3.11) 
we conclude from the above that the union of sequences 

{^•(A,x)} iez , {^(A,x)} iez (3.12) 

for all A of the interior of a also forms a complete system of generalized eigenf unctions of 
the "limiting" selfadjoint operator Jo- (a;), acting in / 2 (Z). The system is known in spec- 
tral theory as the quasi-Bloch generalized eigenf unctions, because in the case of periodic 
coefficients (see e.g. ()3.43J) ) they are well known Floquet-Bloch solutions of correspond- 
ing finite-difference equation. In this context z/(A) is called the quasi-momentum as the 
function of spectral parameter. 

Recall now that if A = {A(x)} x£ Tq-i is a selfadjoint quasiperiodic operator in Z 2 (Z), 
then its Integrated Density of States can be defined as 

MA) = / (^(cc)((A,oo))oo dx, (3.13) 

where £a(x) is the resolution of identity of A(x) (see [T31 I2H] for this and a more general 
case of ergodic operators). 

By using the above definition and ()3.7|) . we find 

k M . ) (d\)=(f V a (X,x)dx)dX, (3.14) 

where kj^^dX) is the measure, corresponding to the non- increasing function &j CT (.)(A) in 

(HOT. 

Another definition of the Integrated Density of States is as follows. Consider the 
restriction A n (x) of A(x) to a finite interval [l,n], imposing certain selfadjoint boundary 
conditions at the endpoints of the interval. The spectrum of A n (x) is a finite set and 
we can introduce its Normalizes Counting Measure of eigenvalues k n as the divided by 
n number of eigenvalues of A n (x) in the interval (A, oo) (cf (jl.4j) \ It is known (see e.g. 
[1211211) that k n converges weakly to (|3.13|) for any x G T^ 1 . 

In the case of operators J a {x), possessing the complete family of quasi Bloch general- 
ized eigenfunctions (j3.5j) - (J3.6J) . it can be shown that &j CT ( ) = v-, i.e. that the Integrated 
Density of States of J a (-) coincides with its quasi-momentum as the function of the spec- 
tral parameter, and we have from ()3.14|) 

k M .)(dX) = v{dX) (3.15) 
= ( / V a (X,x)dx ) dX. 

We found, in particular, a relation between two quantities of asymptotics (|2.8|) . 
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Another important characteristics of quasi periodic (more generally, ergodic) operators 
is the Lyapunov exponent 7a(A), defined as the rate of exponential growth of the Cauchy 
solutions of the corresponding finite difference equation of second order. The Lyapunov 
exponent and the Integrated Density of States are related by the Thouless formula (see 
e.g. [29 , formula (11.82)). In the case of the quasiperiodic Jacobi matrix J a (x) the 
formula is 

7j ct ( )(A) = - / K a {x)dx+ I \og\\ - n\kj^.){d\). (3.16) 

Jjl- 1 Jo 

Since the generalized functions of J a [x) are bounded and do not decay at infinity (see 
(1331) - flUU) or (I3~TTH) - (EHTjn . we have 

7j ct ( )(A) =0, A G a. 

Hence, the l.h.s. of 1)3.16]) is zero if A G a. In view of ()3.14|) the obtained relation is just 
the Euler-Lagrange equation ()2.3|) for the functional ()2.2|) . 

3.2 Orthogonal polynomials with respect to varying weights 

We will present here constructions, similar to those of the previous subsection but for 
orthogonal polynomials with respect to varying weights. To this end it is useful to make 
explicit the amplitude of the potential V in (|1.1U|) . i.e. to replace V by V/g, g > 0. In 
what follows we will keep V fixed and vary g. Thus orthonormal polynomials ffl.l lj) and 
related quantities will depend on g. To control this dependence we will use results of 
papers [T2ll2Tj. 

Note first that if the potential is real analytic, then the minimizer N of 1)2.12]) possesses 
a density p supported on a finite union of finite intervals a [TO] IT1] . 

According to ^S] asymptotics (|2.19|) - (|2.2(J|) are most precise and well behaving if a 
real analytic potential, satisfying (|2.11j) . is regular (see [15] . formulas (1.12) and (1.13)). 
This condition implies, in particular, that the density of the measure N in ()2.13|) is strictly 
positive on the interior of its support cr, and vanishes as a square root at each edge of a. 
Furthermore, following [21j, we say that g is regular for V if V/g is a regular potential. 
If g is regular for V, and 

i 

°go = {J[ai(9o),bi(go)} 
i=i 

is the support of the equilibrium measure N go corresponding to V/go, then there exists an 
open neighborhood of go, consisting of regular values g for V and 

i 

^ = 1>^UG?)] (3-17) 
1=1 

with the same number q of intervals. Besides, ai and hi are real analytic, a\ is strictly 
decreasing and b\ is strictly increasing in g. 

We will also need the following formula, relating N g and v gi minimizing correspond- 
ingly (J2.12|) with V/g instead V and (J2.2|) with a g instead of a jT2] : 

N g =g- 1 [ v g ,dg'. (3.18) 
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The formula was proved in i n a fairly general setting, including piece-wise continuous 
Vs. Its particular cases are given in [2E1HI], where its spectral and asymptotic meaning 
made explicit, related to a kind of "adiabatic" regime in g for corresponding Jacobi matrix 
(ll.20|) (see also the derivation of formula ()4.3|) below). 

Now we can give an analog of constructions of previous subsection, i.e. the "limiting" 
Jacobi matrix with quasiperiodic coefficients. We confine ourselves again to the case of 
rationally independent components of vector (3 of (J2.18|) . a generic case in g. 

Consider the coefficients 7j of the Jacobi matrix (jl.2(Jj) . associated with orthonor- 
mal polynomials {pf^}z>o with varying weight. Introducing explicitly the dependence of 
coefficients on g and writing in view of (jl.lUj) with V/g instead of V 

V V 

n- = l—-, (3.19) 
g gl/n 

we get 

r\ n \g)=rf\gl/n). (3.20) 
Setting here I = rii(x) + k, where now (cf (|3.2J0 

km {rii(x) ft} = x h I = l,..,q-l, (3.21) 

i— >oo 

x = {xi}1~l is a point of T 3-1 , and k is an arbitrary fixed integer, we obtain in view of 
(12.20)) and the continuity of TZy in g and x, and of P in g: 

^r^l(9) = ^v( !! ^^9dn i (x)+k)p(^^ g ]) (3.22) 



n,(,)+*w " ^oo'^V Mx) V'W'-r^^ n^x) 
= Hy (g, ka(g) + x) . 

where 

a(g) = (gP(g))'. (3.23) 
Analogous relation is valid for the diagonal entries of J^ n >: 

lim s nAx)lk(g) = S v (g, ka(g) + x) . 

% — >-oo * v ' 

Now, by using formula f|3.18jl . we find an important relation 

a(g) = a(g), (3.24) 

where a(g) is defined by ()2.6|) with u g instead of v. We conclude from the above that 
limiting coefficients are (cf ()3.3j) ) 



S n ii x )) 

rii(x)+k 

lim &£l k (g) = S v (g, ka(g) + x) , k G Z. 



km r^ x x »(g) = TZ V (g, ka(g) + x) , k e Z, (3.25) 

i— >oo 



As a result we obtain a quasiperiodic Jacobi matrix Jy/ g (x), defined by the coefficients 
(cf(Q) 

IZy (g, ka(g) + x) , Sy (g, ka(g) +x) , k e Z, (3.26) 
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and having the frequencies (ai(g), ...,a q -i(g)), obtained from (|2.5|) - (|2.6j) in which v is 
Ug, hence a is the support a g of N g . Note that Jv/ g (x) is the limit in the sense (|3.21j) in 

j{n 

V/g,W 



the strong operator topology of Z 2 (Z) of matrices Jy) , whose coefficients are defined by 



formulas ()3.1|) with r^_ k and s^_ k instead of r n+k and s n+ k- 

Applying the same limiting argument to asymptotic formula ()2.19|) . and by using 
(|3.21|) and ()3.18|) . we obtain for any fixed fceZ 



1>%U+kW = (2Py(A^,^(^) + x)) 1 / 2 cos(vrn,(x)iV 9 (A) (3.27) 
+ nku g (X) + Q v (\,g, ka(g) + x)j + o(l), n^x) 



oo. 



By using these formulas, the exponential decay of ip^ k (X) outside a g , and the limit 
(j3.21|) . we obtain the complete system of generalized eigenf unctions and the resolution of 
identity £j v/g { x ) of Jv/ g (x), given by (j3.5|) - ()3.6|) (or (j3.12j0 and ()3.7|) . in which subindex 
a is replaced by the subindex V/g. In particular, we have for the diagonal entries of 

£j vh {x)'- 

{ S Jv, 3 {*)) kk ( dX ) = Xa g (X)V v (X,g, ka(g) + x), k e Z, (3.28) 

where Xa is the indicator of a (cf ()3.9|) with j = k). The limit here is the weak limit of 
measures. The support in A of the r.h.s. of this formula is the support a g of the equilibrium 
measure N g . This implies that the spectrum of the quasiperiodic matrix J v / g (x) is o g . 

Note that the spectrum of the "initial" double infinite matrix J„ , defined analogously 
(j3.1j) but via rS fe and s^) k , is M for all n < oo. 

Besides, arguing as in obtaining (|3.14|) . we find for the Integrated Density of States 
measure k Jv/g{ . ) of J V / g (x): 

k Jv/g{ . } (dX) = Vg {dX) (3.29) 

= I / T>v(X,g,x)dx)dX. 
\7t9- 1 / 

Comparing the first equality of this formula with the first equality of (J3.15|) in which a 
is replaced by cr g , we conclude that J Ug (x) and Jv/ g {x) have the same spectrum and the 
same Integrated Density of States. 

The coincidence of spectra of J Ug (x) and Jy/ g {x) implies (see [T7J 1221 EE] ) that each 
of them is an isospectral deformation of another, i.e. that the coefficients f!3.4|) of Ja g {x) 
differ from the coefficients (|3.26|) of Jyj g {x) just by a shift of their argument. This fact 
can also be checked directly, by comparing explicit formulas for both sets of coefficients, 
given [7J El EH] and in ^5] correspondingly, and by using again the trick with infinitesimal 
variation of the amplitude of potential (see Appendix). 

Here is one more link between two classes of polynomials and spectral theory. It 
concerns the Lyapunov exponents of J a {x) and Jy/ g (x) and the potential. It can be 
shown ^2] that the Lyapunov exponents of the both matrices coincide and if j g (X) is 
their common value, then 



V(X)=2 / lgl {X)dg' 1 X G a g . (3.30) 
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3.3 Periodic Jacobi matrices 

Here we consider a class of polynomial potentials V in (11.10)1 for which corresponding 
Jacobi matrices J Ug (x) and J v / g (x) have periodic coefficients. Besides, several quantities 
related to the matrices and orthogonal polynomials can be found explicitly. We follow 

M- 

Let v be a polynomial of degree q with real coefficients and with the leading term z q . 
Assume that there exists g > such that all zeros of the polynomial v 2 — 4g are real and 
simple and set 

VW = (3.31) 
2q 

We will show that in this case coefficients of J ag {x) and Jv/ g (x) are g-periodic and their 
spectrum is 

a g = {A : t> 2 (A) - Ag < 0}. (3.32) 
We show first that the equilibrium measures N g and v g for ()2.12j) and (|2.2)1 are: 

N g (d\) = Pg (X)dX, Pg (X) = l -^\v 2 (X) - Ag\ 1/2 Xa g (X), (3.33) 

and 

u g (dX) = d g (X)dX, d g (X) = ^^ 2 (A) - 4g\~ 1/2 Xa g {X)- (3.34) 
Indeed, it is a matter of direct calculations to find that if p g is given by ()3.33)) . then 

wv/ g (z) : = - \og(z - n)p g (fi)d]i = - u(z)^u 2 (z) - 1 (3.35) 

Ja g 1 

- log (u{z) + y/v?(z) - l) - u 2 (z) - (2q)- 1 \ogg/e. 

where u = v/2^/g, and we use the branch of logarithm with the cut (— oo, 0) and the 
argument n on the upper edge of the cut and the branch of \/u 2 — 1, such that \Ju 2 — 1 = 
u + o(l), u — > oo. 
On the other hand 



$lwv/ g (X + iO) = — / log |A — p\p g (n)dn 

J a a 



is the logarithmic potential of N g (dX) = p g (X)dX. Now, analyzing the values of $ of 
()2.16)1 in this case 

$ 9 (A) = + 2^w v/g {X + tO) (3.36) 

with V from (|3.31j) . we can check directly the validity of (|2.14)1 - (|2.15)1 (with V/g instead 
of V) with the strict inequality in (j2.15j) and ly = (2q)~ 1 \ogg/e. This proves that N g of 
(13.33)) is the minimizer of ()2.12)) with V j g instead of V. 

It can also be proved that v g of ()3.34)) is the minimizer of (|2.2)) with o g of ()3.32)) instead 
of a. We can use either (|3.18j) or the above scheme, computing (cf ()3.35l 



w ag (z) : = -/ \og(z - p)d g (p)dp 

= --log (u(z) + ^u 2 (z) - lj - ^logp, 
13 



and then checking directly (|2.3|) - (|2.4|) . 

We will use, however, another argument to prove (|3.33|) and (|3.34j) . The argument 
is based on a representation, important in the inverse problem for periodic operators of 
second order [2H I2H] ■ 

Let u be a polynomial of degree q with real coefficients and such that all zeros of u 2 — 1 
are real and simple. Then u can be written in the form 

u(z) = cos9(z), (3.37) 

in which 6(z) is the conformal map of the open upper half-plane C + = {z G C : Imz > 0} 
onto the domain 

{6 : qm < m < q 2 ir, %6 > 0} (3.38) 

\ [J {6 : m = In, qi < I < q 2 , < %6 < hi}, 
qi<l<q2, 

Here q± < q 2 are integers, q 2 — qi = q, < hi < oo and 9(oo) = oo. In fact, the r.h.s. 
of ()3.37|) is a polynomial of degree q if and only if — oo < q\ < q 2 < oo, q 2 — q 1 = q |2"4"j . 
Function 6{z) is analytic in C + and continuous in the closed upper half-plane C+. When 
z = A + iO varies from — oo to oo, the limiting value #(A + iO) runs along the boundary 
(the "comb") of the domain ()3.38J) . so that either 9?#(A + iO) varies from {q\+l — l)7r to 
(gi + 1)tt and Q9(\ + iO) = 0, if A varies through the /th "band" [a;, bi], I = 1, q of a g , 
or 3?0(A + iO) = (mod7r) and ^^(A + iO) = k, < k < hi, if A varies through the /th 
"gap" (bi,a t+1 ) of o g . 

By using the terminology of mathematical physics we can say that 6(z)/Trq is an 
analytic continuation of the quasimomentum as a function of energy in the extended 
band scheme. 

We set u = v/2y/g in ()3.37|) . Then the zeros of u 2 — 1 are the band edges — oo < a x < 
bi < ... < a q < b q , 6{b q ) = 0, q\ = —q, q 2 = and #(A + iO) varies from (— q + I — l)n to 
(— q + 1)tt, I = 1, q when A varies from ai to bi in the Ith band. By using (J3.37|) and 
(I3~^|l we can rewrite RZT^ - (l2~TT)|) . as 



$(A) = -ilog,/e+( f °> ifAGK,6,],/ = l, ,g, 
KJ q &y/ \ 0(«i(A)), if A G [b h a l+1 ], I = l,..,q, V ; 



where ai + i = a\, 



. , 2 /sinh2/« \ 4 f K . , 2 , n 
g( K ) := n ) = - / smh tdt > 0, k > 0, 

9 V 2 J q Jo 



and K;(A) varies from and hi > 0, when A varies through the gap (6/,a; + i). This yields 
(|2.14|) - (|2.16|) . thereby a proof that (|3.33j) is the density of the equilibrium measure N g , 
corresponding to the potential ()3.3H) . Moreover, since the inequality in ()2.15|) is strict in 
this case, the corresponding value of g is regular for the potential (|3.31|) . 
It follows from ()3.35|) that 

Ng(X) = -TT-^WV/giX + iO), 

and then ()3.37j) implies that 
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and similarly 

u g {\) = ^ + (A), (3.41) 

where 9+(X) = 3R6 1 ( A + iO), and 9(z) is defined in ()3.37|) . In view of the above properties 
of this function, we have 

N g {a l+1 ) = v g {a l+1 ) = 1—t, 1 = 1,..., q-l, (3.42) 

and then (J2~HJ) and (12~T%D imply 

a, = A = — l = l,-,g-l. (3.43) 

Hence, the coefficients (13.4)1 of the matrix J ag (x) and the coefficients (|3.2fi|) of the matrix 
Jy/ g (x) are g-periodic in this case. Moreover, we need not to consider in this case the 
whole torus T q ~ 1 as the set of values of x in (|3.2|) and (j3.21|) . but just the set of vertices 
of the regular g-polygon. This is similar to a standard procedure of the theory of almost 
periodic functions, where the corresponding set is the closure of all limiting points of 
sequences (13. 2 j) or (|3.21j) . hence depends on arithmetic properties of the frequency vector 
a or (3. 

Consider simple cases of potentials (|3.31|) . The case q = 1 corresponds to v(X) = —A 
and yields 

<r g = [-2y/g,2y/g\, 

P 9 (A) = ^(4^-A 2 ) 1 /V 9 (A), 



d 9 {\) = \Ag-X 2 )- 1 l\ ag (X). 



The first density corresponds to the well known semicircle law by Wigner for the Gaussian 
Unitary Ensemble The role of polynomials pf 1 ^ play hi(X^n/2g)(n/2g) 1 ^ A , where 

{hi}i>o are the orthonormal Hermite polynomials. The second density is the Density of 
States of the Jacobi matrix with constant coefficients 77 = ^Jg, si = 0, I 6 Z. The matrix 
plays here the role of both limiting matrices J ag and Jy/g ■ 

The case q = 2 corresponds to v(X) = A 2 + Vo, vq < —2^fg, and yields 

a 9 = [- b (9),-a(9)} U [a{g),b(g)\, 
a{g) = (\v \ - 2^~gfl\ b{g) = (\v \ + 2^) 1/2 , 

P9(A) = 2^^ 2 - A2)(A2 - a2) ) 1/2x - (A) ' 

d 9 (A) = ^((6 2 -A 2 )(A 2 -a 2 ))- 1/2 XCT9 (A). 

Asymptotics of corresponding orthogonal polynomials were considered in jS]. Matrices 
J ag and Jy/g are both of period 2 and their Density of States is given above. 

For a general two interval case, where the corresponding matrices are quasiperiodic and 
their coefficients can be expressed via the Jacobi elliptic functions see (21 El E2] (ordinary 
polynomials) and |Hj (polynomials with varying weights). 
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The fact that in the case of ordinary polynomials the limiting finite band Jacobi matrix 
is periodic if its spectrum is the inverse image of a polynomial map (see (|3.32Jl ) is known 
(see e.g. j^U EE! an d references therein). It is of interest that the same property holds 
also for polynomials with varying weights and that the corresponding potential (J3.31|) is 
also polynomial and can be explicitly related to the map. 

We mention one more link of asymptotics of orthonormal polynomials and periodic 
Jacobi matrices ^JJ that concerns the Hill discriminant (or the Lyapunov function) of 
Ja g (x) and Jy/ g (x) and the polynomial v of (j3.31|) . Recall that the Hill discriminant A (A) 
is defined as 1/2 of the trace of the monodromy (transfer) matrix of corresponding finite- 
difference equation of second order with periodic coefficients and plays an important role 
in spectral theory (see e.g. [221 EH])- It can be shown [TTj that both matrices have the 
same Hill discriminant A g and that 

A, = v{X)/2y/g. 

We discussed above the case, where the polynomial v in ()3.31|) is such that all zeros of 
v 2 — Ag are real and simple. Admitting non-simple (but still real) zeros, we include the 
case, where two adjacent bands touch one another or a band is going to appear inside a 
gap. 



4 Eigenvalue distribution of random matrices 

Here we discuss briefly certain aspects of eigenvalue distributions of ensembles (JUT} 
(I1.3J1 . related to the above results, in particular to the matrix Jyi x y 

4.1 Expectation of linear statistics 

According to (TUl [Tj3] we have for any bounded and continuous <p 



lim E{N n [<p]} = [ <p{\)N{d\), 

J a 



where N is the minimizer of (|2.12jl . Combining this with (|3.18jl . (|3.28jl . and I3.29|) . we 
obtain for the r.h.s. of this formula 

<p(\)d\ / v g (d\)dg = dg (tp(Jy/ g (x))) 00 dx. 
Jo Jo Jit- 1 

On the other hand we can always write (|1.6|) as 

N n [<p) = n- lr Tr<p(M n ), 
and we obtain a kind of "functional correspondence" 

lim E{n- 1 Tr^(M n )} = f dg [ { V {J v/g {x))) m dx, (4.1) 

reminiscent to that for ergodic operators, see ^H|, Theorem 9.6, and [2E], Theorem 4.4. 

Here is a heuristic argument, explaining the above formula. According to ()1.14|) the 
l.h.s. of the formula includes the orthonormal functions ^ of ()1.12j) for I = 0, ...,n — 1. 
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Indicating explicitly the dependence of these functions on g and using the relation (cf 

ri n \\,g)=4 l \\,gl/n), (4.2) 
we obtain from (|2.19|) that the leading contribution to p n as n — > oo is 



n— 1 

-Vzv(A,Z/rM/3(Vn)) 



1=0 



Assuming that T>y(X, g,x), and /3(g) are continuous in g, we can say that the summand 
in this formula is "slow varying" ml/n and "fast varying" in I. This observation results 
in the limiting formula for the density p(A) of the measure N: 

p(X) = [ dg [ V v (X,g,x)dx. (4.3) 
Jo Jji- 1 

Using now (I3~23J) and (JTTHjl . we obtain flUpJ . 

4.2 Covariance of linear statistics of eigenvalues 

By using (|1.21j) we write (jl.l5|) as 

Cav{N n [<p 1 ],N n [<p 3 ]} = ^ f f A ^ A ^ C n (X 1 ,X 2 )dX 1 dX 2 , 

where Aip/AX is defined in p.lfijl . 

C n (A 1; A 2 ) = (r^f (e^^AOei^^Aa) -e£_i(Ai)e2_i(A 2 
and 

e£(A)=W"' ) (A)^ ) (A) 

is the density of the resolution of identity {£j(^)(dX)}'^ iri=0 of matrix (jl.20j) (cf (|3.7Jl ). 
Thus, assuming that ipx,2 are bounded and of the class C 1 and passing to a subsequence 
{rii(x)}i>i that satisfies (|3.21|) . we obtain in view of (|3.22|) 

lim n i (x) 2 Cov{N ni(x) [(p 1 },N niix) [ip 2 \} = Kl(x-a) (4.4) 

m(xj— >oo 

^T^f ( e °'°^ 1 ^ e - 1 -iCV) ~ e o ,-i(Ai)e -i(A 2 ))rfAidA 2 , 

where ej t k(X), j,k G Z is the density of the (jk)th entry of the resolution of identity 
{£j v (x)(dX)}i t m € z of the limiting Jacobi matrix Jy(x), determined by (|3.26|) . 

This result seems rather unusual from the point of view of traditional probability 
concepts. Indeed, the covariance of linear eigenvalue statistics (11.6)1 is of the order rT 2 
rather than of the order rT x as in the case of independent identically distributed random 
variables, or in a more complex and close to our context case of the Schrodinger operator 
with random potential. In the latter case it is a matter of routine spectral theory argument 
to show that, say for 

<p s (\) = (\-z)- 1 , SszjLO, (4.5) 
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where 

N n [<p g ] =n- 1 Tr(H n -z)- 1 

and H n is the discrete Schrodinger operator on the interval [l,n] with a random i.i.d. 
potential, then 

lim nCovln^Tr^ - z x )-\ n^Tr^ - z 2 )- 1 } = C{ Zl ,z 2 ), 

where C(zi, z 2 ) is analytic for Qzi >2 7^ and is not identical zero if the second moment of 
the potential exists. 

On the other hand, it follows from (|4.4|) that if M n is a random n x n random matrix, 
given by ffTTTl) - (Ol) . then 

lim n J 2 (x)Cov{n i (x)~ 1 Tr(M n - zj) -1 , n i (a;)- 1 Tr(M n - z 2 )- x 

n;(a;)— >oo 

_ ^2 / ^ AG 0i o(x) AG-x-i(x) fAG ^(xy 2 ' 



Az Az V Az 



where 



{( J V(x) - Z\) - \-Jv(x) - Z 2 ) ), 



Az z\ — z 2 J ' 

Hence the covariance of N n [(p z ] = n -1 Tr(M„ — z)~ l is of the order 0(n~ 2 ). 

This indicates that the Central Limit Theorem, if any, should be valid not for n 1 ' 2 N n \(p z 
as in the case of i.i.d. random variables, but for nN n [ip z ], i.e. for the sum 



(nh 

1=1 



without a n-dependent factor in front. 

This was indeed shown in for the single interval case q = 1 and for a rather broad 
class of test functions. However, as we have seen above, the case q — 1 is exceptional, 
since it is only in this case asymptotic formulas (j2.19|) and ()2.20j) do not oscillate in n 
because of the absence of the argument n/3 in corresponding coefficients of the formulas. 
Hence, for q > 2 the limiting normal law for J\f n [(p x ], if it exists, could be different for 
subsequences in (13. 2j) having different limits in TP 7-1 , because its variance depends on 

lim Var{J\f ni ( x )[ip z ]} 



where 



1Zy(x - a)^g 0fl (z,x)g- 1 - 1 (z,x) - \g ,-i(z,x)\ 2 \ 



9jk{z,x) 



\X-z\ 2 

However, as is shown in [2H], the situation with the Central Limit Theorem for linear 
eigenvalue statistics of random matrices (jl.2j) - (jl.3|) is more subtle. Namely, the above 
scheme of a family of the Gaussian limiting law with the x-dependent variance for various 
subsequences in (|3.21j) proves to be valid in the case, where the matrix Jy{x) is periodic. 
In a generic case of quasiperiodic Jy{x) the limiting laws of subsequences {Af ni (x) M}i>i 
exist but are not Gaussian. 



18 



A Appendix 



Here we verify directly that the coefficients of the Jacobi matrices J ag (x) and Jv/ g {x) 
coincide up to a shift in x. We will consider again the generic case, where the frequencies 
()2.6j) and (j2.18|) are rationally independent, hence x varies over the whole T q ~ l . Besides, 
we consider only the off-diagonal entries of J ag (x) of ()3.4|) and Jy/ g (x) of ()3.26|) (note that 
the diagonal entries are zero if V is even). 

Recall that the coincidence follows also from general results on the inverse problem of 
spectral analysis for "finite-band" potentials, known as the algebro-geometric approach 
(see e.g. |T7J 123 EH] ) • Indeed, since the spectra J ag (x) and Jy/ g (x) coincide (see Section 
3.2), J ag (x) is a isospectral deformation of Jv/ g (x) and vice versa, hence, by the inverse 
problem, one of them can be obtained by a shift in x of another [TT| l3~H]). 

We begin a direct proof of this assertion by recalling necessary results of spectral 
theory of finite band Jacobi matrices and related facts of complex analysis on Riemann 
surfaces (see e.g. [T71 I22J EE|)- 

Given the set o of (j2.1|) . denote T the two sheeted (hyperelliptic) Riemann surface, 
defined by the equation 

q 

w 2 = R(z), R{z) = \{{z - ai)(z - &,), 
i=i 

i.e. obtained by pasting together two copies of the complex plane along the union of the 
"gaps" (bi, a 2 ), (65-1, a q ), (b q , a x ) of a, the last gap goes through the infinity point. Let 
idp be the normalized differential of the third kind with simple poles of residues ±1, at the 
infinity points P± on each sheet of T, and let U = (U\, U q -i) be the vector of 6-periods 
of dp: 

Ui = ±- [ dp, Z = l,...,g-1, (A.l) 
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2tt 

where {biYiZi are the so-called 6-cycles on T. 
On the other hand, the integral 

-p 

idp, P ,P G r 



Po 

with a properly chosen initial point P can be identified with the complex Green function 
G(z) of C \ a with the pole at infinity (see e.g. jJJ). The real part g(z) = $t.G(z) is 
uniquely determined by the requirements to vanish for its limiting values on a and to be 
harmonic in C \ o for g(z) — log \z\. It follows then that if v is the unique minimizer of 
(12 .2|) . hence solves the corresponding Euler-Lagrange equation (J2.3j) . then 

g(z) = / log \z - n\v(dfi) - l a /2. 

J s 

This and (jA.l|) imply (|2.6|) . where aci — = N([ai+i, I = 1, •••?<? — 1 is the 

harmonic measure at infinity of the (/ + l)th "band" [aj+i, of a. 

Denote 9 : T 9_1 — > C the Riemann ^-function, associated with V. Then according to 
[101 ES] the leading coefficient 7„ of the polynomial p n , where {pi}i>o are orthonormal 
polynomials on o with respect to weighs, satisfying (J2.7j) . is for n —>■ 00: 

9{na + tt(oo) + d a ) 



9{na — u{oo) + d a 



+ o(l 



(A.2) 
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Here l a is defined in (|2.3|) . 

= I lu 



with the integral taken along a path on the first sheet and uo = (coi, to q --i) is the 
canonical basis of the differential of the first kind on T, A a and d a do not depend on n 
but depend on cr, the weight, and the points d,...,(q-i of T that are the poles of the 
corresponding Baker- Akhiezer function [TJj. In the case, where ai = mi/q with positive 
integers mi, ...,m q , hence with a (/-periodic J a (x) (see e.g. Section 4), ([, ■■■,C q -i are the 
eigenvalues of the Dirichlet problem on the period for the corresponding finite-difference 
equation, distributed in a fixed way over the edges of the gaps. These are in fact the 
parameters, indexing representatives of the isospectral family. Another characterization 
of ([, C_i is given in [7j, Theorem W2. 
Asymptotic formula (jA.2|) and the relation 

r n = lnlln+i, (A. 3) 

expressing the off-diagonal entries of a Jacobi matrix via the leading coefficients of asso- 
ciated orthonormal polynomials, lead to the relation 

r 2 _ c -i a 9(( n + ty» - Moo) + d a )6(na + u(oo) + d a ) | 
9{{n + \)oc + u(oo) + d a )9{na — u{oo) + d a ) 

Replacing here n by n + k, where k is an arbitrary fixed integer (in fact k = o(n)), and 
passing to the limit (|3.2j) . we obtain for the function lZ a of ()3.4j) : 

7^ ( x ) _ e -ia @( x + a - u (°°) + d a )9(x + u(oo) + d a ) ^ 
9(x + a + u(oo) + d a )9(x — u(oo) + d a ) 

By using the formula 

a + 2m(oo) = 0, (A.6) 
that follows from the Riemann bilinear relations (see e.g ^Zj, Section 6), we can write 

n ff {x) = TZ{x + x a ), (A.7) 

where 

. , , 9(x + a)9(x — a) , k , 

n{x)=e ekk ~ (A ' 8) 

and 

x (j = —u(oo) + dfj . (A. 9) 

Consider now the orthonormal polynomials {p\ n ^}i>o with respect to varying weights (jl.lOj) 
- (jl.lip . Then we have for the leading coefficient of according to j!5j . formula (1.63): 



9(nf3 + u(oo) + d v ) 

9{n(3 - u(oo) + d v ) + °^ \ 



n -> oo, (A.10) 



where Zy is defined in (I2.14|) - (|2.1bjl . (3 is defined in (|2.18|) . u{oo) is the same as in (jA.2j) . 
and Ay and dy do not depend on n but depend on V and the points Cg-u that are 

zeros of a certain analytic function on C \ cr (see [T3], formulas (1.26) - (1.27), (1.30)). 
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In view of the relations 

W _ (») / (») c A ii\ 

(cf ()A.3p ) and (|3.25|) we need the coefficients 7^7 fc , n — > oo, 6 Z fixed (see ()3.22|) ) 
in order to find the entries of Jy(x). We will find them by using the same trick as in 
obtaining ()3.22j) . According to the trick the passage from n to n + k, n — > oo, A; = o(n) 
can be carried out by passage from the super-index n to n + k, which is equivalent to 
the infinitesimal change g —>■ g + gk/n in the inverse amplitude of the potential. Thus, 
replacing V by V/g in (jA.10|) . using the above trick and (|3.23j) - (|3.24|) . we obtain (cf 



6{nf3 + ka + u(oo) + dy) ^ 
#(n/3 + ka — u(oo) + dy) 



n — > oo, 



Now, comparing the Euler-Lagrange equations (EHl) - (l2~TT)J) for (l2~T2l and (Q - (Q 
for O, we find the relation (cf (l3~23]) - lCT23j) ): 



k s = {gW/g)'- 

This and ()A.11|) yield for the asymptotics of the off-diagonal entries of matrix j( n ) (jl.2(J|) . 

\l>0- 

2 / , , \ 2 



associated with {p[^}i 



-2 



_^ 6*(n/3 + {k + l)a — -u(oo) + dy)Q{nf3 + ka + w(oo) + (iy) 
~~ 6 0(n/3 + (k + l)a + u(oo) + d v )9(n/3 + ka - u(oo) + d v ) + °^ 

Passing here to the limit (J3.21j) . we obtain for the function IZy of (|3.25|) . determining the 
off-diagonal entries of the limiting matrix Jy(x): 

lZy(ka + x) := TZy(l, la{l) + x) = 1Z{ka + x + xy), 



i.e., 



K v {x) = K{x + Xy) (A.12) 
where TZ is defined in (TOl) . and (cf ijOjl ) 

Xy = — «(oo) + C?y/ 9 . (A. 13) 

Comparing ()A.7|) and ()A.12|) we conclude that lZ a and 7£y differ by a shift of argument. 
The assertion, formulated at the beginning of the Appendix is proved. 
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